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ABSTRACT 

We study the anisotropy of Ultra-High Energy Cosmic Ray (UHECR) events collected by the 
Telescope Array (TA) detector in the first 40 months of operation. Following earlier studies, we 
examine event sets with energy thresholds of 10 EeV, 40 EeV, and 57 EeV. We find that the 
distributions of the events in right ascension and declination are compatible with an isotropic 
distribution in all three sets. We then compare with previously reported clustering of the UHECR 
events at small angular scales. No significant clustering is found in the TA data. We then check the 
events with E > 57 EeV for correlations with nearby active galactic nuclei. No significant correlation 
is found. Finally, we examine all three sets for correlations with the large-scale structure of the 
Universe. We find that the two higher-energy sets are compatible with both an isotropic distribution 
and the hypothesis that UHECR sources follow the matter distribution of the Universe (the LSS 
hypothesis), while the event set with E > 10 EeV is compatible with isotropy and is not compatible 
with the LSS hypothesis at 95% CL unless large deflection angles are also assumed. We show that 
accounting for UHECR deflections in a realistic model of the Galactic magnetic field can make this 
set compatible with the LSS hypothesis. 
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1. INTRODUCTION 

One of the keys to understanding the nature of the 
Ultra-High Energy Cosmic Rays (UHECRs) is their dis- 
tribution over the sky. This distribution depends on the 
UHECR sources, as well on the UHECR mass composi- 
tion and large-scale magnetic fields, both Galactic and 
extragalactic. Despite significant effort, none of these 
issues is well understood at present. 

Observation of the cutoff in the highest - energy 
part of the cosmic ra y spectrum ([Abbasi et al.l [2008a; 
lAbraham et al.ll2008"B ) suggests that the UHECR propa- 
gation at high energies is limited by the interaction with 
the cosmic backgroun d radiation (the Greisen-Zatsepin- 
Kuzm in (GZK) effect (|G rciscn 1966; Zatscpin & Kuzmin 
1966)). One therefore expects that the closest sources of 
UHECRs are situated within the GZK volume of the size 

100 Mpc. At these scales the matter distribution in the 
Universe is inhomogeneous, and so must be the distribu- 
tion of the UHECR sources. If propagation of UHECRs 
at these distances is quasi-rectilinear (whether or not this 
is the case depends on both their composition and the 
magnetic fields), one generally expects the UHECR flux 
to be anisotropic, showing variations at large angular 
scales and possibly point sources. 

If UHECR primary particles are protons, as suggested 
by the composition measurements performed by the High 
Resolution Fly's E ye (HiRes) and the Telescope Arra y 
(TA) experiments (jAbbasi et al.ll2010bt iTamedal [2010h . 
the UHECR propagation is, in fact, expected to be 
quasi-rectilinear. Wit h the existing estimates of the 
Galactic magnetic f ield (Han ct al. 2006; Sun et al.ll2008l : 
iPshirkov et alJ 20111) and bounds on the extragalactic 
ones (jKronbere 119941 ) , the deflections of protons should 
be relatively small. For instance, a random extragalac- 
tic field of magnitude 1 nG and correlation length of 
~ 1 Mpc would deflect a proton of energy 10 20 eV by 
about 2° over a distance of 50 Mpc, while the Galactic 
field would produce deflections of order 2 — 4° depend- 
ing on the direction. In this case a sizeable anisotropy 
may be expected regardless of the density of the UHECR 
sources down to energies as low as 10 19 eV. 

On the contrary, if the composition at highest ener- 
gies is heavy or predominantly heavy, as the results of 
the P ierre Auger Observatory (PAO) ([Abraham et alJ 
12010( 1 seem to indicate, the quasi-rectilinear propaga- 
tion is not expected for the bulk of UHECRs. Some 
anisotropy at large angles may still arise if the ex- 
tragalactic fields are sufficiently small and the den- 
sity of sources is such that only a few nearby ones 
contribute to the observed flux, but the small-scale 
anis otropy would b e supp r essed (for recent ana lyses see, 
e.g.. iGiacinti et al.l ([21)1(1 ; iTakami et alJ ([20121 )'). Thus, 
the study of the UHECR anisotropy may shed light 
on both the ma ss composition and the density of the 
UHECR sources dDubovskv et alllgOOOHYoshiguchi et alJ 
120031 [2001 IKachelriess fc Semikozll2005l )~ 

Numerous attempts at detection of the UHECR 
anisotropy have been made previously. Early studies 
indicated clustering of the UHECR events at small an- 



gular scales ([Havashida et al.|[l996t iTinvakov fc Tkachevi 
120011) . On the basis of small-scale correlations, differ- 
ent classes of pu tative sources o f UHECR were sug- 
gested (see, e.g.. IGorbunov et al.l (120041) : lAbbasi et alJ 
(2006); lAbraham et al.1 ([20071 l2008aD ). More recently, 
the Pierre Auger Observatory has claimed correlations 
of UHE CRs with the nearby Active Galactic Nuclei 
(AGN) ([Abraham et al.1 120071 l2008al ) which were not 
confirmed by observ ations in the Northern hemisphere 
([Abbasi et al.ll2008bf ) . 

At larger angular scales, evidence for correla- 
tion s with the s uperg alac tic plane was claimed 
bv IStanev et al.1 (ll995l). IGlushkov! ([200l . and 
IGlushkov fc PravdinI (I2001D but not confirmed by 
other authors (IHavashida et al~l 119961: IKewlev et al.1 
[19961: iBuTetaLl I1999D . Also. IKashti fc Waxmanl ([20081 ) 
found the anisotropy in the PAO data which was not 
confirmed by the HiR es data in the Northern hemisphere 
([Abbasi et al.ll2010al) . 

In this paper we present the anisotropy analysis of 
UHECR observed by the Surface Detector (SD) of the 
Telescope Array in the first 40 months of its operation. 
TA is a hybrid UHECR detector located in the Northern 
hemisphere in Utah, USA (39°17'48" N, 112°54'31" W) 
which has been fully operational since March 2008. It 
consists of 507 scintillator detectors c overing the area of 
approx imately 700 km 2 (for details see lAbu -Zavvad ~et al.1 
(|2012bf )). The atmosphere over the surface array is 
viewed b y 38 fluorescence tele scopes arranged in 3 sta- 
tions (see lTokuno et al.l ([2012T )). The surface detector of 
TA is the largest in the Northern hemisphere. 

In this paper we focus on testing previous observations. 
Namely, we consider the clustering of the UHECR events 
at small angular scales (as would be produced by bright 
point sources) , possible correlation with nearby AGN and 
correlation of the TA events with the large-scale struc- 
ture (LSS) of the Universe. Following previous analy- 
ses, we consider three a priori chosen energy thresholds: 
10 EeV, 40 EeV, and 57 EeV. It should be noted that dif- 
ferent experiments may have different energy scales due 
to different systematic errors in the energy determina- 
tion, which may affect the selection of the events. When 
referring to the results of other experiments, we assume 
the energy scales as reported by these experiments. In 
statistical tests which require a pre-defined confidence 
level we set the latter to 95%. 

The paper is organized as follows. In Section [2] we 
describe the data sets used. In Section [3] we examine 
the event sets for a presence of small-scale clustering by 
studying the UHECR auto-correlation function. In Sec- 
tion [4] we consider correlations of UHECR events with 
nearby AGNs. Section [5] describes our search for correla- 
tions of the UHECR events with the large-scale structure 
of the Universe. In Section [6] we summarize the results 
and present conclusions. 
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Among the existing TA data sets (SD data, Fluores- 
cence Detector (FD) data in mono and stereo modes, 
and hybrid detector data) the SD data set has by far the 
largest number of events. The present analysis is based 
on the data collected in the period 2008.05.11-2011.09.15 
(40 months) of operation by the TA surface detector ar- 
ray. Cutting events with zenith angle > 45°, the SD data 
set contains 988 events with energies > 10 EeV, 57 events 
with E > 40 EeV, and 25 events with E > 57 EeV. This 
is the largest UHECR set to date in the Northern hemi- 
sphere. 

The angular resolution of TA events with E > 10 EeV, 
is approximately 1.5°. This follows from the comparison 
of the thrown and reconstructed arrival directions of sim- 
ulated data sets, and is supported by the direct compar- 
ison between the SD and FD arrival directions of hybrid 
events. The energy resolution of the TA surface detector 
at E > 10 EeV is better than 20% (|Abu-Zavvad et al.l 
I2012al) . 

The exposure of the TA surface array is calculated by 
the Monte-Carlo technique with full simulation of the 
detector, which will be described elsewhere. As follows 
from the Monte-Carlo simulations, the acceptance of the 
TA surface detector for E > 10 EeV and zenith angle 
cut of 45° is close to the geometrical one. For reasons of 
computational efficiency, in the present analysis aimed at 
anisotropy at relatively small angular scales we use the 
geometrical acceptance to generate random event sets. 

Figure [1] shows the comparison between the distribu- 
tions in declination (left column) and right ascension 
(right column) of the events simulated with the geomet- 
rical exposure (red line) and the data (blue data points) 
at the energy thresholds of 10 EeV, 40 EeV, and 57 EeV 
(top, middle, and bottom rows, respectively). The com- 
patibility of expected and observed distributions in all 
6 cases was checked by the Kolmogorov-Smirnov (KS) 
test. The lowest KS probability was p = 0.13 for the 
distribution in the right ascension at E > 57 EeV. Thus, 
all three sets are compatible with a uniform distribution. 

3. AUTOCORRELATION FUNCTION 

The AGASA experiment reported clustering of 
UHE CR events with E > 40 EeV at the angular scale of 
2.5° (fHavashid a et al.lll996() . Here we repeat this analy- 
sis using the TA data set. 

The procedure is as follows: for a given angular separa- 
tion, 5, we count the number of pairs of observed events 
that are separated by an angular distance less than 5, 
thus obtaining the data count. We then generate a large 
number (typically, 10 5 ) of Monte-Carlo (MC) event sets 
each having the same number of events as the real data 
set. The simulated sets are generated with a uniform dis- 
tribution according to the TA exposure. In each MC set 
we count pairs of events in the same way as in the data, 
which gives the MC count for that set. We then calculate 
the average MC count for all of the MC sets. This repre- 
sents the expected number of pairs for the angular scale 
5, assuming a uniform cosmic ray distribution. For each 
value of 5 we then determine the fraction of simulated 
sets where the number of pairs is greater than or equal 
to the number of pairs in the data. This gives the p- 
value, P(5), that describes how likely the excess of pairs, 
if found in the data, is to occur as a result of a fluctuation 
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Fig. 1. — Comparison between the data (blue points) and the 
sets of 10 4 events simulated with the geometrical exposure (red his- 
togram) at energies 10 EeV, 40 EeV, and 57 EeV (top, middle, and 
bottom rows, respectively). Plots show the distribution of events in 
declination (left column) and right ascension (right column). The 
compatibility of the two distributions by the Kolmogorov-Smirnov 
test is given as Pks in the upper left corner of each plot. 

in a random set. Small values of P(S), thus indicate a 
departure from uniformity at the corresponding angular 
scale. 

We first perform a blind test of the AGASA claim. 
Fixing the energy threshold to 40 EeV and the separation 
angle to 5 = 2.5° we find pairs while 1.5 pairs are 
expected in the case of a uniform distribution. Therefore, 
there is no excess of small-scale clusters in the TA data. 

We next extend the analysis to all angular scales. No 
significant excess is found. The results are illustrated 
in Figure [2] for angles from to 40° and three energy 
thresholds of 10 EeV, 40 EeV, and 57 EeV as specified 
on the plots. For each energy threshold, the left panel 
shows the number of pairs with the angular separations S 
binned in 2° bins (data points) . The shaded region repre- 
sents the average number of pairs expected in the case of 
the uniform distribution. Both the data and the uniform 
expectation are normalized bin-by-bin to the area of the 
bin, so that in the case of a uniform full-sky exposure 
the expectation would be flat. The overall normalization 
is set in such a way that the expectation in the first bin 
equals one. 

The right panels of Figure [2] show the dependence of 
the p- value, P(S), on the separation angle, S, for the 
corresponding energy. Note that P(S) is a cumulative 
quantity since it takes into account all the pairs sepa- 
rated by angles from to S. For this reason a small, but 
coherent over several bins, excess at angles from 10 to 
20 degrees on the lower left panel of Figure [5] produces a 
more significant feature in the corresponding -P(<5), lower 
right panel of Figure [5J This feature corresponds to the 
group of events visible on the sky map (see the lower 
panel of Figure [7]). 
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Fig. 2. — Autocorrelations in the TA data sets at E > 10 EeV, 
E > 40 EeV, and E > 57 EeV (top, middle, and bottom rows, 
respectively). Left panels: the number of pairs with angular sepa- 
rations <5 normalized to the area of the angular bin (data points), 
compared to the expectation for the uniform distribution (shaded 
histogram). The errors are 1-sigma Poisson errors. Right panels: 
probability, P(S), that the excess of pairs with the angular separa- 
tion less than 8 occurs as a fluctuation in a uniform distribution. 
Small P(S) indicates a departure from isotropy. 

When accessing the significance of departures from 
isotropy on the basis of P(S) represented in Figure [21 
one should take into account the fact that the angular 
scale of the excess is not known in advance. Thus, there 
is a stat istical penalty for choosing this scale a posteri- 
ori (see iTinvakov fc Tkachevl ([2004) for a detailed dis- 
cussion). Taking this penalty into account, none of the 
three examined data sets shows a significant deviation 
from an isotropic distribution. 

Interestingly, although close clusters in the high-energy 
TA event set are absent, one of the TA events falls within 
1.7° of a h igh energy event o bserved by the Auger Ob- 
servatory (jAbreu et al.l I2010h . Both events have E > 
10 20 EeV. The center of the doublet has the Galactic 
coordinates I = 36°, b = —4.3°. 

4. CORRELATION WITH ACTIVE GALACTIC NUCLEI 

The Auger collaboration h as reported a correlation 
(|Abraham et al.l 120071 l2008al) between UHECRs with 
E > 57 EeV and the nearby (redshift z < 0.018 or, equiv- 
alently, distance d < 75 Mpc) Active Galactic Nuclei 
AGNs) from the Veron-Cet ty & Veron (VCV) catalog 
Veron-Cettv fc Veronl 12006). The greatest correlation 
was observed at the angle of 3.1°. In the control data set, 
the number of correlating events was 9 out of 13, which 
corresponds to about 69% of events. The Auger col- 
laboration has recently updated the analysis and found 
that a smaller fraction of the UHECR events correlates 
wit h the same set of A GNs in the latest UHECR data 
set (jAbreu et al.l l20l"oh than in the original one. Out of 
55 events with E > 55 EeV, 21 were found to correlate 
with AGNs, which corresponds to a fraction of correlat- 
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Fig. 3. — Hammer projection of the TA cosmic ray events with 
E > 57 EeV and nearby AGNs in the Galactic coordinates. Corre- 
lating and non-correlating events are shown by filled red and empty 
blue circles, respectively. AGNs are represented by black dots. The 
dashed line shows the boundary of the TA exposure. 

ing events equal to 38%. In this section we test the TA 
data for correlations with AGN. 

The set of 472 nearby AGNs used by lAbraham et al.l 
(2007) contains 7 objects listed at zero redshift, all in 
the field of view of TA. Of these 7 objects, two are stars, 
one is a quasar with unknown redshift, one is a Seyfert 2 
galaxy, two are spiral galaxies (including the Andromeda 
galaxy) and one is a dwarf spheroidal galaxy. We exclude 
these objects from the analysis, which leaves 465 objects 
in the AGN catalog. 

The TA exposure is peaked in the Northern hemi- 
sphere, so that the AGNs visible to TA are largely dif- 
ferent from those visible to Auger, though there is some 
overlap. The distribution of nearby AGNs over the sky is 
not uniform because of the large scale structure (see Sec- 
tion [3] for more detail) and because the VCV catalog is 
not complete: due to observational bias it tends to con- 
tain more objects in the Northern hemisphere. For this 
reason, a larger fraction of events is expected to correlate 
with AGNs in the TA data under the assumption that 
AGNs are sources of the observed UHECRs. Taking into 
account the distribution of nearby AGNs over the sky 
and assuming equal AGN luminosities in UHECR, we 
estimated the correlating fraction will be ~ 73% for TA 
on the basis of the original PAO claim, and ~ 43% on 
the basis of the updated analysis by PAO. 

The sky map of TA events with E > 57 EeV and 
nearby AGNs from the VCV catalog is represented in 
Figure [3] in Galactic coordinates. The cosmic rays are 
shown by filled red (correlating events) and empty blue 
circles (non-correlating events). AGNs are shown by 
black dots. 

Figure [4] shows the number of TA events correlating 
with AGNs as a function of the total number of events 
with E > 57 EeV ordered according to arrival time. 
The black dashed line represents the expected number 
of random coincidences in case of a uniform distribu- 
tion calculated via Monte-Carlo simulation. The blue 
line shows the expected number of correlating events as 
derived from the original PAO claim. Shaded regions 
represent 68% and 95% CL deviations from this expec- 
tation _cah2ukj : ted_j2ythfi_jnaximum likelihood method of 
Ref. (|Gorbunovet al.l [20061. As is seen from Figure [U 
present TA data are compatible with both isotropic dis- 



Search for Anisotropy with TA 



5 




5 10 15 20 25 



W,„, 

Fig. 4. — The number of TA events with E > 57 EeV correlating 
with VCV AGNs as a function of the total number of events. The 
expectation according to the original PAO claim is represented by 
the blue line together with the 1- and 2-sigma significance bands. 
The black dashed line shows the expected number of random co- 
incidences. 

tribution and the AGN hypothesis. 

In the full TA SD data set, there are 11 correlating 
events out of 25 total, while the expected number of ran- 
dom coincidences for this total number of events is 5.9. 
Making use of the binomial distribution with the proba- 
bility of a single event to correlate pi SO — 0.24, one finds 
that such an excess has probability of ~ 2% to occur by 
chance with isotropic distribution of arrival directions. 

5. CORRELATION WITH LSS 

Even though the sources of UHECRs are not known, 
their distribution in space at large scales must follow 
that of the ordinary matter. The latter is anisotropic at 
scales below ~ 100 Mpc forming the Large-Scale Struc- 
ture (LSS) of the Universe that consists of galaxy clus- 
ters, filaments and voids. If UHECRs are not strongly 
deflected on their way to Earth, their distribution over 
the sky should correlate with the nearby structures, with 
over-densities corresponding to close clusters and under- 
densities corresponding to voids. 

The amplitude of anisotropy depends on the UHECR 
propagation length (the larger is the propagation length, 
the smaller contributions of the local structures and, 
therefore, the anisotropy) and on the UHECR deflec- 
tions. In this section the propagation of UHECR is cal- 
culated assuming they are protons. However, it should 
be noted that regardless of whether the UHECR compo- 
sition is heavy or light, their propagation length changes 
with energy roughly in the same way and becomes of or- 
der several tens of megaparsecs as the energy approaches 
10 20 eV. Thus, the most important parameter that de- 
termines the amplitude of the anisotropy at a given en- 
ergy is the typical deflection angle which we denote as 
9 (which is, of course, very different for heavy and light 
composition). 

The goal of this analysis is to determine which values 
of 6 are compatible with the space distribution of the TA 
events. In principle, this can be done at all energies. To 
minimize statistical penalties, we limit our analysis to 
the energy thresholds of 10 EeV, 40 EeV, and 57 EeV. 

5.1. Statistical Method 

To test the compatibility between the observed 
UHECR distribution over the sky and that expected 



under the LSS hypothesis (that is, the hypothe- 
sis that UHECR sources trace matter distribution in 
the Universe), we employ the method developed by 
IKoers fc Tinvakovl ()2009bf ) and used previousl y in the 
analysis of the HiRes data (Abbasi ct al. 2010a). In this 
method, one first computes the UHECR flux distribution 
expected under the LSS hypothesis and then compares 
it to the observed one by the flux sampling test. 

The matter distribution in the nearby Universe may 
be inferred from the complete galaxy catalogs containing 
the redshift information. In this work we use the 2MASS 
Galaxy Redshift Catalog (XSCzQ that is derived from 
the 2MASS Extended Source Catalog (XSC), with red- 
shifts that have either been spectroscopically measured 
(for most of the objects) or derived from the 2MASS 
photometric measurements. This catalog provides the 
most accurate information about 3D galaxy distribution 
to date. 

For the flux calculations, we use the flux-limited sub- 
sample of galaxies with apparent magnitude m < 12.5. 
For fainter objects, the completeness of the catalog 
degrades progressively, while their inclusion does not 
change the results considerably. We exclude objects 
closer than 5 Mpc in order to avoid breaking of the sta- 
tistical description (if such objects are assumed to be 
sources of UHECR, they have to be treated individu- 
ally). We also cut out galaxies at distances further than 
250 Mpc replacing their combined contribution by a uni- 
form flux normalized in such a way that it provides the 
correct fraction of events as calculated in the approx- 
imation of a uniform source distribution. The quanti- 
tative justification of thes e procedures can be found in 
(|Koers fc Tinva kov 20093). The resulting catalog con- 
tains 106 218 galaxies, which is sufficient to accurately 
describe the flux distribution at angular scales down to 
~ 2°. The UHECR flux distribution is reconstructed 
from this fl ux-limited catalog b y the weighting method 
proposed by Lvnden-Belll (11971f) and ad apted to flux cal- 
culations by IKoers fe Tinvakovl (|2009af) . 

The XSCz catalog loses completeness in the band of 
roughly ±10° around the Galactic plane and especially 
around the Galactic center. The size of this region is not 
much larger than a typical deflection of a proton even 
at 57 EeV, so this gap may be bridged without loss of 
accuracy. Away from the Galactic center at \l\ > 60° 
where only a fraction of the galaxies (the dimmer part) 
is missing in the catalog, we apply a I— and b— dependent 
weight correction to the remaining galaxies so as to com- 
pensate for the missing ones. In the region close to the 
Galactic center, \l\ < 60°, we extrapolate the flux den- 
sity from the adjacent regions in a straightforward man- 
ner. The latter is not an accurate procedure; however, 
the Galactic center region overlaps with the TA exposure 
only slightly, and this inaccuracy is not important for our 
results as can be checked by excluding this region from 
the analysis. 

When propagating the UHECR primary particles from 
a source to the Earth, we assume them to be pro- 
tons and take full account of the attenuation processes. 
The injection index at the source is taken to be 2.4, 
which is compatible with the UHECR spectrum observed 

1 We are grateful to T. Jarrett for providing us with the prelim- 
inary version of this catalog. 
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Fig. 5. — Sky map of expected flux at E > 57 EeV (Galactic 
coordinates). The smearing angle is 6° . Letters indicate the nearby 
structures as follows: C: Centaurus supercluster (60 Mpc); Co: 
Coma cluster (90 Mpc); E: Eridanus cluster (30 Mpc); F: Fornax 
cluster (20 Mpc); Hy: Hydra supercluster (50 Mpc); N: Norma 
supercluster (65 Mpc); PI: Pavo-Indus supercluster (70 Mpc); PP: 
Perseus-Pisces supercluster (70 Mpc); UM: Ursa Major (20 Mpc); 
V: Virgo cluster (20 Mpc). 

by HiRes and TA (jAbu-Zavvad et alj[2012ah assuming 
proton composition and the source evolution parameter 
m = 4 (|Gelmini et al.l[2007t ). We also assume that the 
effects of both the Galactic and extragalactic magnetic 
fields can be approximated by a single parameter, the 
Gaussian smearing angle 9. We consider 9 a free param- 
eter and vary it in the range 2 — 20°. In general, the 
deflections of UHECR in magnetic fields contain both 
random and regular parts, the latter being due to the 
regular component of the Galactic magnetic field. The 
regular deflections are not Gaussian. However, the sta- 
tistical test we use here is not sensitive to the coher- 
ent character of deflections provided they do not exceed 
10 — 20° as set by the typical size of the flux variations 
due to local structures (cf. Figure [5]). So, for the most 
part of the analysis we will use the Gaussian smearing 
to represent all the deflections without making the dis- 
tinction between the regular and random ones. Later, 
in Section 15.41 we will discuss, in the case of the lowest 
energy set and the largest deflections, the effect of explic- 
itly accounting for the regular component of the Galactic 
magnetic field. 

To calculate the expected flux, we assume that 
UHECR sources follow the space distribution of galax- 
ies. The simplest way to realize this assumption in 
practice is to assign each galaxy an equal luminosity in 
UHECR's. This is a good approximation if the density 
of the UHECR sources is sufficiently high (so that many 
sources are present in local structures contributing to the 
anisotropy). The contribution of each galaxy to the total 
flux is then calculated taking into account the distance 
of the source and the corresponding flux attenuation. 
Individual contributions are smeared with the Gaussian 
width 9, so that the flux at a given point of the sky is a 
sum of contributions of all the galaxies within the angu- 
lar distance of order 9. Further deta ils on the flux cal- 
culation can be found in References (iKoers fe Tinvakovl 
I2009bt iKoers fe Tinvakovll2009at lAbbasi et al.ll2010aD . 

Figure [5] shows the flux map calculated by the above 
procedure for an energy threshold of 57 EeV and smear- 
ing angle 9 = 6°, not yet modulated with the TA expo- 
sure. Darker regions correspond to higher flux. A band 
of each color integrates to 1/5 of the total flux. One 
can identify the nearby structures which are marked by 
letters on the plot as explained in the caption. 



The next step is to compare the calculated flux dis- 
tribution to the actual distribution of the TA events 
and determine whether they are statistically compati- 
bl e. In this work we use the f lux sampling test proposed 
bv IKoers fe Tinvakovl (|2009b| ). The starting point is the 
map of the expected flux calculated as explained above. 
One reads off the flux values at positions of the data 
events. This gives a set of numbers which we refer to as 
the "data set" . One may say that the cosmic ray events 
sample the flux map in a particular way that depends 
on their space distribution. One then generates a large 
number of Monte-Carlo events which are distributed ac- 
cording to the expected flux and reads off the flux values 
at their positions. This gives the set of flux values which 
we refer to as the "MC set". If the angular distribu- 
tion of the data and MC events is the same, so must be 
the distributions of the flux values in the data and MC 
sets. These two distributions may be compared by the 
parameter-free Kolmogorov-Smirnov (KS) test. 

The result of the KS test is the p-value which shows 
whether the data and MC flux sets are drawn from the 
same parent distribution. If this p-value is low, the two 
distributions of flux values are different and, therefore, 
the angular distributions of data and MC sets arc differ- 
ent. 

5.2. Estimate of Statistical Power of the Flux Sampling 

Test 

An important characteristic of a statistical test is its 
ability to discriminate between two hypotheses, or the 
statistical power. For the case at hand, the statistical 
power is the probability to rule out the LSS hypothesis 
at 95% CL if the cosmic ray distribution is isotropic. The 
closer the statistical power to one, the more sensitive is 
the test. Knowing the statistical power provides an a 
priori idea of what kind of sensitivity can be reached 
with the given number of events. 

In general, the statistical power increases with smaller 
smearing angles since this improves the contrast in the 
flux map. For the same reason, the statistical power 
increases with energy (the UHECR propagation length 
becomes shorter and the relative contribution of the local 
structures is therefore enhanced). Also, the statistical 
power increases with the number of events. 

We have calculated the statistical power of the flux 
sampling test in case of TA for the three energy thresh- 
olds of 10 EeV, 40 EeV, and 57 EeV, and smearing angles 
varying from 2° to 14°. We have found that for the ac- 
tual number of events in the TA data set, the statistical 
power is below 50% for smearing angles 9 > 9°, 9 > 3° 
and 9 > 4° for the above three energy thresholds, respec- 
tively. 

The case E > 57 EeV is shown in Figure [6] The vari- 
ous curves in the plot correspond to different number of 
events (note that the actual number of events in the TA 
data set with E > 57 EeV is 25). The gray region rep- 
resents the expected range of deflections in the Galactic 
magnetic field in the case of protons. 

5.3. Results 

First, we check the compatibility of the TA event sets 
with the isotropic distribution. To this end we generate 
an isotropic flux map modulated with the TA exposure. 
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Fig. 6. — The statistical power of the flux sampling test at 
E > 57 EeV as a function of the smearing angle. Different curves 
correspond to different number of events, as indicated on the plot. 
The actual number of events in the TA data set with E > 57 EeV 
is 25. The gray region shows the expected range of deflections in 
the Galactic magnetic field in the case of protons. 

This map is independent of energy and smearing angle. 
We then test the compatibility of the TA event set for 
E > 10 EeV, E > 40 EeV, and E > 57 EeV with this 
map. The flux sampling test gives the p-values 0.5, 0.9 
and 0.6 for the three data sets, respectively. Thus, at all 
three energy thresholds the data appear to be compatible 
with an isotropic distribution. 

Next, we examine the compatibility of the TA event 
sets with the LSS hypothesis. Figure [7] shows the 
skymaps of the expected flux at energy thresholds of 
10 EeV, 40 EeV, and 57 EeV (top to bottom) and the 
smearing angle of 6°. The white dots represent the ar- 
rival directions of the TA events. The bands are drawn 
in the same way as in Figure [5l i.e., each band integrates 
to 1/5 of the total flux. This means, in particular, that 
if the LSS model were true each band would contain 1/5 
of the total number of events in average. Note that the 
configuration of the bands changes with energy because 
of the energy dependence of the propagation length. 

The results of the flux sampling tests are presented in 
Figure |U The p-values are shown as a function of the 
smearing angle at energy thresholds of 10 EeV, 40 EeV, 
and 57 EeV. Each point represents the p-value obtained 
by the flux sampling test at the corresponding energy 
threshold and smearing angle. 

As one can see from the plots, for E > 40 EeV and 
E > 57 EeV the data are compatible with the structure 
hypothesis at the 95% C.L. The decrease of the p-values 
slightly below the 95% C.L. in the case E > 57 EeV can- 
not be assigned a real significance in view of the penalty 
factors for trials (e.g., three energy thresholds). 

Although large smearing angles do not have a straight- 
forward physical interpretation in view of the Gaussian 
approximation used, we have investigated the behavior 
of the p-values corresponding to the case E > 57 EeV 
for larger smearing angles and found that it fluctuates 
around p ~ 0.05 for angles as large as 9 ~ 50° and then 
goes to p ~ 1. Such behavior may arise because the flux 
map for E > 57 EeV remains anisotropic even for very 
large smearing angles. 

At the energy threshold of E > 10 EeV the situation 
is somewhat different. The data are incompatible with 
the structure model up to angles of order 20°. In view 
of the large deflections in magnetic fields at low ener- 
gies, such behavior is expected. One should be careful, 




Fig. 7. — The skymaps of the expected flux at energy thresholds 
of 10 EeV, 40 EeV, and 57 EeV (from top to bottom) in Galactic 
coordinates with the TA events superimposed (white dots). The 
smearing angle is 6°. 

however, with the interpretation of this result. First, 
Figure [8] does not include the penalty for the number of 
trials. Second, at E > 10 EeV the uncertainties in the 
flux calculation due to the choice of the model parame- 
ters (in particular, the injection index and the evolution 
parameter) are the largest. Finally, if the smearing an- 
gle is attributed to deflections in the magnetic fields, the 
dominant contribution is likely to come from the regular 
component of the Galactic magnetic field, as discussed 
in the next section. Such large and regular deflections 
require a more accurate modeling, which we attempt in 
the next section. 

5.4. Accounting for the Galactic Magnetic Field 

The deviation from the structure model at E > 10 EeV 
and small smearing angles is an indication that magnetic 
field deflections play an important role in the distribu- 
tion of the UHECR arrival directions. In general, several 
contributions to the deflections are expected. First, there 
are deflections produced by intergalactic magnetic fields. 
These fields are known quite poorly. They are usually 
thought to obey the upper boun d of B <, 10~ 9 G with 
a correlation length, I <J 1 Mpc (Kronberg 1994). With 
these parameters, a proton of energy 10 EeV coming from 
50 Mpc would be deflected by ~ 20°. However, there are 
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Fig. 8. — The results of the statistical test for the compatibil- 
ity between the data and the LSS hypothesis. The p-values (red 
points) arc shown as a function of the smearing angle 9. Low p- 
values indicate incompatibility with the LSS model. The horizontal 
line shows a confidence level of 95%. The three panels correspond 
to energy thresholds of 10 EeV, 40 EeV, and 57 EeV from top to 
bottom, as indicated on the plots. 

indications that the extragalactic mag netic fields may be 
several orders of magnitude smaller (jDolag et al.l l2005h 
than the upper bound. 

Second, UHECRs are deflected in the regular compo- 
nent of the Galactic Magnetic Field (GMF). The regular 
GMF is known much better than extragalactic fields. It 
can be inferred, e.g., from the Faraday rotation measures 
of Galactic and extragalactic radio sources. According 
to recent studies, a typical deflection of a 10 EeV proton 
would be 20 - 40° (Ps hirkov etaLl 120111) . This is com- 
parable or larger than the deflection in the extragalactic 
field. 

Finally, the Galactic field has a random component. 
Although the amplitude of this component is a few times 
larger than the regular one, its contribution into the de- 
flections is subdominant (or at most comparable) to that 
of the regular component (|Tinvakov fc Tkachcv [2005ft 
due to its random character. 



From this discussion it is clear that, most likely, the 
regular part of the magnetic field provides the dominant 
contribution into the UHECR deflections. At low en- 
ergies when the magnitude of the deflections becomes 
large, Gaussian smearing is not a good approximation 
for such deflections. They have to be taken into account 
explicitly. 

In order to see whether the deflections in the regular 
GMF can be a reason for the discrepancy between the 
data and the LSS model we have repeated the analysis 
of Section [5751 including the regular GMF. The presence 
of the regular magnetic field is taken into account by 
modifying the expected flux distribution. The smearing 
angle remains a free parameter; it accounts for random 
deflections in the extragalactic fields and in the random 
component of the GMF. The statistical test itself remains 
unchanged. 

We adopt the recent GMF model by I Pshirkov et al.l 
(2011). This model has been obtained by fitting the 
GMF model parameters to the latest catalog of the Fara- 
day rotation measures of extragalactic sources. In addi- 
tion to the disk field, this model also contains a toroidal 
halo field. 

Although the fits to the Faraday rotation measures 
constrain the parameters of the GMF, some combina- 
tions of these parameters are constrained rather poorly. 
In particular, the magnitude of the halo field is degener- 
ate with the halo height above the Galactic disk: making 
the halo field stronger and simultaneously higher above 
the disk does not strongly affect the rotation measures. 
Thus, there remains some freedom in the choice of the 
GMF parameters. The question is whether this freedom 
can be used to bring the arrival directions of UHECR 
into accord with the LSS hypothesis without contradict- 
ing the Faraday rotation data. 

We have found that the compatibility with the LSS 
model cannot be reached without the halo field. When 
the halo is included, the compatibility with the LSS 
model is possible, although the required halo field is 
rather strong (but still compatible with the data on the 
Faraday rotation measures) . 

An example of the flux map with the GMF included is 
shown in the upper panel of Figure [§] The flux distribu- 
tion is calculated for the case E > 10 EeV and smearing 
angle of 6°. The magnetic field parameters are as follows: 
the magnitude of the halo 4 /J.G and the thickness of the 
halo is 1.5 kpc. Note that after the inclusion of GMF the 
Virgo region has moved away from the TA field of view, 
and the expected flux distribution has become closer to 
the uniform one. 

The results of the flux sampling test of the LSS model 
with the regular GMF included are shown in the lower 
panel of Figure HO Black squares represent the p-values 
in the case of the GMF with the parameters described 
above. For comparison, red circles show the p-values 
in the absence of GMF (the same as the upper panel 
of Figure [8]) , while green triangles represent the case of 
GMF with the disk component only. One can see that the 
regular GMF can produce deflections that make the data 
for E > 10 EeV compatible with the LSS model for all 
but the smallest smearing angles. Thus, the discrepancy 
between the LSS hypothesis and the TA data with E > 
10 EeV can, in principle, be explained by the deflections 
in the regular GMF. 
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Fig. 9. — Upper panel: The sky map of the expected flux for 
E > 10 EeV and smearing angle 6° taking into account the GMF 
(Galactic coordinates). The parameters of GMF are as follows: 
the magnitude of the halo is 4 fiG and the thickness of the halo 
is 1.5 kpc. Note the absence of overdensity in the direction of 
the Virgo cluster. Lower panel: The result of the statistical test 
of the compatibility between the TA event set with E > 10 EeV 
and the LSS hypothesis for different models of GMF: no magnetic 
field (circles), disk component only (triangles), both disk and halo 
components (squares). The horizontal line shows the confidence 
level of 95%. Low p-values indicate incompatibility. 

6. SUMMARY AND CONCLUSIONS 

In this paper we present a search for anisotropy in the 
TA data collected over the period of about 40 months, 
which is the largest UHECR data set to date in the 
Northern hemisphere. The main focus of this paper is 
on checking the existing claims: small-scale clustering, 
correlation with nearby AGNs, and correlation with the 
LSS. 

The results are summarized as follows. 

• The TA data show no clustering of the UHECR 
events at small scales, neither at the angular scale 
of 2.5° in the set with E > 40 EeV as reported by 
the AG AS A experiment, nor at any angular scale 
from to 40 degrees in the data sets with E > 
10 EeV, E > 40 EeV, and E > 57 EeV. There 
is a hint of grouping of events at angular scales 
of 20 — 30° at the highest energies, however the 
statistical significance of this feature is insufficient 
for a definite conclusion. 

• There is no statistically significant correlation of 
the TA data with E > 57 EeV with the posi- 
tions of nearby AGNs from the VCV catalog us- 
ing the parameters reported by the PAO (angular 
scale of 3.1° and redshift cut in the VCV cata- 
log z < 0.018). Out of 25 observed events with 



E > 57 EeV, 11 have been found to correlate with 
positions of nearby AGNs, while 5.9 are expected 
in average from random coincidences (chance prob- 
ability of 2%). 

• The TA event sets with E > 10 EeV, E > 40 EeV, 
and E > 57 EeV appear compatible with a uniform 
distribution according to the flux sampling test. 
The sets with E > 40 EeV and E > 57 EeV are 
also compatible, at 95% CL, with a model which 
assumes that sources follow the large-scale struc- 
ture of the Universe (LSS model). The set with 
E > 10 EeV is not compatible, at 95% CL, with the 
LSS model unless the deflections of these UHECR 
exceeds 20°. 

• The set with E > 10 EeV can be made compatible 
with the LSS model, at smearing angles larger than 
~ 3°, by including the effect of the regular compo- 
nent of the Galactic magnetic field and assuming a 
realistic model for the latter. The smearing angle 
in this case represents the deflections in the random 
Galactic and extragalactic fields. 

From the analysis presented, one concludes that there 
is no apparent deviation from isotropy in the present TA 
data. At high energies, this may be merely due to an 
insufficient number of events. However, if this tendency 
persists at several times larger statistics, it will be diffi- 
cult to reconcile with the proton composition of UHECR 
regardless of the source nature: if the sources within the 
GZK volume are numerous, they must follow the (inho- 
mogeneous) matter distribution. If the source density 
is small so that there are only a few within the GZK 
volume, this very fact will produce anisotropy. 

At lower energies the deflections are expected to be 
large even for protons, which makes the distribution of 
the events more isotropic. However, the number of events 
is much larger as well, so that even small deviations 
from isotropy may become detectable as the statistics 
increases. The fact that for E > 10 EeV the distribu- 
tion of the events is not compatible with the LSS model 
without assuming a large (il 20°) smearing angle may 
indicate that we are observing the first manifestation of 
UHECR deflections in the Galactic magnetic fields. The 
possibility to reconcile the observed UHECR distribution 
with the LSS model by correcting for the deflections in 
the realistic model of GMF is in accord with this inter- 
pretation. 
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8. APPENDIX: LIST OF EVENTS WITH E > 57 EEV 

In this Appendix we present the list of events with 
energy E > 57 EeV and zenith angle 9 < 45° that have 
been recorded by the surface detector of the Telescope 
Array from May 11, 2008 to September 15, 2011. During 
this period, 25 such events were observed. Table Q] shows 
the arrival date and time of these events, the zenith angle 
9, energy in units of EeV, and Galactic coordinates I and 
b in degrees. The angular resolution of these events is 
~ 1.5°, while the energy resolution is better than 20%. 
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